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Properties in magnetic ordered states of graphene nanoribbons with zigzag shaped edges are in- 
vestigated by applying mean-field approximation to the Hubbard model with on-site repulsion U. We 
observe that magnetic moments and critical temperature show anomalous power-law dependences as 
a function of U; the actual values of the power are determined by only the width of ribbons. Such 
' singular behaviours are found to be due to localized nature of the electronic states close to Fermi 

Jo ' energy. 
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as follows. 1 7) The zigzag GNR has a metallic band structure irrespective of the width N in the sense 
that the energy gap does not appear at the Fermi energy in the absence of doping, E = 0. However, 



1. Introduction 

Graphene-based materials with nano-meter sizes have been attracting much attention due to their 
possibilities as new potential devices in application as well as novel stages for emergence of exotic 
phenomena in fundamental science. Especially, the graphene nanoribbon with zigzag shaped edges, 
Tj - ' which is abbreviated to zigzag GNR in the following, is known to have fascinating peculiar properties 

o" 
in 
p 

unlike usual metals, the asymptotic form of the energy dispersion near E = is written as E oc 
±\k — ir/a\ N ,^ and the Fermi velocity vanishes for N > 2 where a and N express the lattice spacing 
and the width of the ribbon (see Fig.l), respectively. Such characteristic properties are due to the fact 

•i-H . 

^ . that the one-particle states near E = are well localized around zigzag edges, i.e., the states close 

to E = are so called edge states. The zigzag edges and the localized states around them have been 
observed by scanning tunneling microscopy and spectroscopy. ^ 

Magnetic properties of the zigzag GNR have been investigated by applying mean-field approxima- 
tion to the Hubbard model with on-site repulsion U 2 ' 14 ~ 17 ) and by the first principles calculation. 18 ~ 22) 
It has been found that the large spontaneous magnetic moments appear at the zigzag edges, which is 
originated from the edge states. The magnetic moments align ferromagnetically at each edge but with 
the opposite direction between the edges. In addition, the magnetic order appears under the infinites- 
imal on-site repulsion; the conclusion is different from the case for graphene sheets where the finite 
amount of U is necessary for emergence of the antiferromagnetic state. The spin excitations 14 ' 22) and 
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the interedge superexchange interaction 17) have been calculated based on the magnetic structure in- 
troduced above. The treatments beyond the mean-field approximation, in which quantum fluctuation 
is fully taken into account, have been carried out. 23 ~ 25) The ground state is found to be Mott insulator 
with charge gap, and the field theoretical approach demonstrates that the Heisenberg model on the 
zigzag GNR expressing the spin excitation belongs to the same universality class as spin 1/2 square 
ladders 23 ^ (gapped for even number legs, gapless for odd number of legs). The result is confirmed by 
numerical calculation. 25 ^ 

In ordered states seen in electron systems, usually, the order parameter shows exponential depen- 
dence as a function of the coupling constant; this fact is due to the finite density of states (DOS) at 
the Fermi energy. On the other hand, the DOS of the zigzag GNR show divergence at Fermi energy. 
Therefore, in the magnetic ordered state of the zigzag GNR, unusual dependence of the spontaneous 
magnetic moments as a function of U is expected. In the present work, we study properties of the 
magnetic ordered states in the zigzag GNR by applying the mean-field approximation to the Hubbard 
model. We focus on interaction dependence of spontaneous magnetic moments and critical tempera- 
ture. It is found that these quantities show anomalous power-law dependences and the actual values 
of the power are determined by only the width N of the ribbons. We clarify that the unusual band 
structure close to Fermi energy gives rise to such singular behaviours. 

2. Model 

The graphene nanoribbons with zigzag shaped edges we study are illustrated in a schematic way 
in Fig. 1. Figs. 1 (a) and (b) show the zigzag GNRs with the width N being even and odd, respectively. 




Fig. 1. Schematic illustration of zigzag GNRs consisting of TV legs with TV = even (a) and N — odd (b). 
Here, the rectangle written by the dash dotted line shows the unit cell and a is the lattice spacing. The filled 
(open) circles express the A (B) sublattices. The two types of slices in the unit cell are denoted as a and /3. 

In the respective figures, the rectangle written by the dash dotted line and a express the unit cell and 
the lattice spacing, respectively. The filled (open) circles indicate the A (B) sublattices. In the unit cell, 
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there are two types of carbon slices; those are denoted as a and /3. 
The Hamiltonian we consider reads % = %\ + Hint, where 

N 



1,(7 1=1 
JV-1 



*E E {4 B ,a( I H+mA I ) + ^ c -} > 



i> i=l 
JV 



(1) 



(2) 



^int E E { n i^v(-O n »>w(-0 + «is,<7(-f)nis,-a(^)} , 

1,(7 1=1 

where i is the hopping between the nearest neighbor carbon atoms. Here c\ x a (I) is the creation 
operator of an electron with the spin <t = ±(<t = +/— expresses f / 4- s pi n state) at the iX site in 
the I-th cell, and riix,a(I) = c\ x a {I)cix,a{I) (* = 1, • • • , and X = A, 73 unless explicitly noted 
in the following). The mean-field approximation is applied to %- mt as 



U 



N 



Hint -> H^ t F ~ CJm ^) C L,a( / ) C ^^( / ) 

7,(7 1=1 

+ (n iS - °miB)c\ Bc (I)CiB,a{I)}, (3) 

where nix = Y2a( n iX,(r(I)) an d m iX = &{nix,a{I)) are the charge and spin order parameters 
with (• • • ) being the thermal average. We should note that the magnetic solution as well as the para- 
magnetic one in the neutral system has uia = = 1 due to the particle-hole symmetry. Therefore, 
in the following, we neglect the terms including the charge order parameters for simplicity since those 
renormalize the chemical potential. As a result, the Hamiltonian under the mean-field approximation 
is written as follows, 

U M¥ = X>t(A:,<7) {M*0 + O)} (4) 

k,a 

where &(k,a) = (c\ Aa (k), c\ Ba {k), ■ ■ ■ , c* NA a (k), c* NB (r (k)). Here h k (k) and hf^(a) are the 
2N x 2N matrices, 
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hfot( a ) = x dia § ( m iA, miB, • • • , mNA,mNB) , (6) 
where jk = 2cos(fca/2). Here, the Fourier transformations 

CiA,M) = 4p E e ifc( - +( - 1)ia/4) ^,,(fc), (7) 



are introduced with Nl being the total number of the unit cell in the system and xj = la. We note 
that h\{k) becomes the real and symmetric matrix owing to the choice of the Fourier transformation, 
eqs. (7) and (8). The 2N order parameters miA and miB are determined self-consistently by 

2N , ^2 

1 \ - 



1 f r i 1 

rEE ff fw]^, f(4 j) m, (9) 

k,a j=l 

2N 2 

^ = ^EE ff {Fw]J /(4%)), do) 



where /(-E 1 ) is the Fermi function defined by f(E) = l/{ex.p(E/T) + 1} with T being the tempera- 
ture. Here E&\k) is an eigenvalue of the matrix h^(k) + /i^ t F (cr) and the corresponding eigenvector 
is expressed by (k), the l-th (I = 1, 2, ■ ■ ■ , 2A) element of which is a real number and written as 
vi^ (k) . The energy per an atom e(U) is given by 

2N N 

e(U) = 2NJV EE^W/^W) + gjv E + . (11) 

fc,CT J=l 1=1 

3. Results and Discussions 

We can obtain the two kinds of self-consistent solutions; one expresses the antiferromagnetic (AF) 
state which satisfies = — m^-i+i)B an< l the other is the ferromagnetic (F) state with miA = 
m(jv-i+i)B- 17 ' 19 ' 21) The magnetic moment at each site and the energy difference Ae(U) = e(U)—e(0) 
in the both states are shown in Figs. 2 (a) and (b), respectively, for the N = 3 system at the absolute 
zero temperature. In each figure, the quantities in the AF state and those in the F state are expressed by 
the solid and dotted curves, respectively. Though difference between TO2A and m^A (miB and m2&) 
in the AF state is not clearly seen in Fig. 2 (a), there exists the significant difference between the both 
quantities. As we expected, the AF state is more stable than the F state. 21) Then, we concentrate on 
the AF state unless noted. Note that the energy difference per an atom between the AF state and the 
F state €af(U) — £f(U) becomes smaller with increasing the width N. This fact seems to indicate 
that difference in the energy is mainly originated from that in the local spin structure, for example, the 
spin configuration at the bond in the center of the ribbon, i.e., nA-nB bond for N = 2n — 1 and nB- 
(n + 1)A for N = 2n. Qualitative discrepancies between the N =odd zigzag GNR and the N =even 
one have been observed in transport properties 26-33 ) and persistent currents. 34) Such discrepancies are 
not found in the magnetic ordered states obtained by the mean-field approximation. 
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U/t 



Fig. 2. Self-consistent solutions of the AF state and the F state for the N — 3 zigzag GNR; magnetic moments 
at each atom in the unit cell, and m^, (a) and energy difference Ae(J7) = e(U) — e(0) (b) where 
e(C7) is defined in eq. (11). In each figure, the solid and dotted curves express the quantities in the AF state 
and those for the F state, respectively. The inset in (b) shows the energy difference between the AF state 
and the F state, £af{U) — cf(U), for several choices of N. 

We investigate in detail the magnetic moment at the zigzag edges, which is far largest than the 
others, and then considered to dominate the magnetic properties as far as U/t <C 1. The magnetic 
moment at the 1A site miA of the AF state is shown as a function of U/t close to U/t = in Fig. 3 
for N = 2, 3, ■ ■ ■ , 10. We can see that the magnetic moment miA behaves as miA oc 
i.e., the quantity shows power-law dependences and the actual value of the power is determined by the 
width N. The result seems to be anomalous because the order parameters of the usual ordered states 
realized in the electronic systems are known to show exponential dependences as a function of the 
coupling constant for weak coupling limit. Actually, in the case of N = 1 which is nothing but the 
usual one-dimensional system, txi\a = —mis — 16(t/U) exp(— 2-irt/U) as long as U/t ^ 1. 

Next we investigate the critical temperature T c of the AF state, which is also known to show 
exponential dependence as a function of the coupling constant in the usual ordered state. Fig. 4 shows 
the critical temperature as a function of U/t close to U/t = 0. Here, the power-law dependence is 
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Fig. 3. Magnetic moment at the 1A site, m\x, of the AF state as a function of U/t for the several choices of 
N. The solid lines express the fitting by miA cx (f//t) 1 /( JV_1 ). 




Fig. 4. Critical temperature T c of the AF state as a function of U /t for the several choices of N. The solid 
lines express the fitting by T c /t oc (f//t) JV /( JV - 1 ). 

also seen in the critical temperature as T c oc 

([//t)^/(JV-i)_ The actual value of the power in T c is 
different from that in mix- However, it can be well understood by the simple dimension analysis, 
T c ~ U x mix, with considering that the AF transition is dominated by the edges states. 

Here, we explore the critical temperature from divergence of the magnetic susceptibility. The 
instability toward the AF state and that toward the F state will be investigated. In the presence of the 
external magnetic field Hix, which does not depend on an index of the unit cell I but does depend on 
the location in the unit cell, the linear response theory results in the magnetic moments, miA and mts 
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as 



N 
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where Hf% = HiX + U rriix/2 is the effective magnetic field at the iX cite. The quantities x^x i'X' = 
lim ? _^o li m cj n -i-o Xix vx' (?) l0J n) (X, X' = A, B) is the susceptibility of the non-interacting system: 
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with Ey\k) (j = 1, • • • , 2iV) being the eigenvalue of fok(&) and v^\k) being the eigenvector corre- 
sponding to it. Namely, the element ofv^\k) expresses the amplitude of the j-th eigenf unction in the 
transverse direction in the non-interacting case. 

In the F state, the magnetic moments satisfy the configuration = rri(N-i+i)B> whereas = 
— TZ(AT-j+i)B i s realized in the AF state. In order to lead to such a magnetic structure, the effective 
external field should satisfy Hf% = H^_ i+1)B ( Hf% = ~H^_ i+1)B ) for the F (AF) state. Based 
on the above consideration, eqs. (12) or (13) are rewritten as follows, 



rfiA = Xf/afHa, 



-l 



Xf/af - [ 1 - -^Xf/af ) Xf/AF' 



(18) 
(19) 



where ttla = (ttiia, fniAi ■ ■ ■ ^naY, Ha = (Hia, H2A, • • • ,HnaY and 1 is the N x N unit 
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(20) 



where the upper and lower sign correspond to [x F ] w anc ^ [^af] a> > respectively. The transition tem- 
perature is determined by det(l — {U/2)x F / AF ) = 0. 

We focus on the transition temperature for weak repulsion, in which case it is numerically obtained 
that the transition temperature is proportional to U N ^ N ^^ . In this case, we take account of only the 
two kinds of eigenstates whose eigenvalues are around E = 0, i.e., E = E^ N+1 ^ = —E^ ~ 
Ajq\k — ir/a\ N with An being a positive constant. Eigenf unctions of such states are known to be 
well localized around the zigzag edges, which are explicitly shown in Fig. 5. At the sufficiently low 
temperature, owing to the localized nature of the eigenf unctions, the matrix elements of Xp/AF can ^ e 
neglected except [x^/^Jn, which are calculated as follows, 

N L ^ de 



r~0i z °J 



e=E, 



(N+l) 



(*) 



[x 



1 



2 ^/(e)-/(-e) 



N L 
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k N+ V(k) Jo I 2e 



(21) 



D N (e). (22) 



Here, Djv(e), which is the density of states close to e = 0, is obtained as Djv(e) ~ Cn/\z\ 1 4) 
from the asymptotic behaviour of the energy dispersion, Eq N+1 ^ = —Eq N ^ = A^\k—TT/a\ N . In deriv- 



ing eqs. (21) and (22), the amplitude at the zigzag edges are approximated as Vq ' (k) = ±1/ a/2 and 
v^^ (k) = ±1/ \/2 (J = N, N+l), and the others are discarded. In addition, the sign of two kinds 



2N 



of amplitudes are assigned according to the result in Fig. 5, e.g., in the odd N case, 



v o 
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1/V2, 



v o 



(k) 



2N 



(k) 



2N 



-l/y/2 from Fig. 5 



= -1/V2, [4 N \k)\ i = -1/^/2, and 
(a). The transition temperature T c f/af> which is determined by 1 — (U/2)[x F / AF ]n = 0, is obtained 
as follows; 



T C , F = (C N X FtN U) N « N -V , 
T c ,af = (CnXaf,nU) n ^ n 1 - ) 



where Xp /af,n is an constant depending on the width N, 

1 1 
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(23) 
(24) 

(25) 
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Fig. 5. (Color online) Amplitudes of the eigenfunctions in the transverse direction close to E/t = for 
N = 5 (a) and TV = 6 (b) where ka/ir = 0.98 is used. In the respective figure, the upper and lower graph 
express the wavefunction of the conduction band and of the valence band, respectively. The slices a and f3 
in the unit cell are defined in Fig. 1. 



f°° 1 1 1 sinhx 

^af,tv = / «x . . (26) 

Jo e x + 1 e x + 1 x v l >" x 

Thus the transition temperatures of the both magnetic states are proportional to U N /( N ~ 1 ^; the result 

for T Cj af is identical with that obtained by the numerical calculation. Note that T c af > T C; p because 

of Xaf,n > ^f.as which corresponds to the fact that the antiferromagnetic state is more stable than 

the ferromagnetic state. 

Finally, we discuss the case where the hopping integral at the 1st leg and that of the iV-fh one are 
modified as c x t and d x t, respectively. In this case, the matrix elements of h^(k) in eq. (5) are 
changed as [h k (k)]i t2 = [h k (k)] 2 ,i = -dr/k and [h k (k)] 2 N-i,2N = [hk(k)] 2 N,2N-l = -dtj k - Even 
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upon such a change, the power of the energy dispersion chose to E = is not changed; i.e., the energy 
dispersion is given by Eq N+1 ^ = —Eq N ^ ~ An(c, c')| k — ir/a\ N as is shown in Fig. 6. The power-law 




Fig. 6. Energy dispersion EQ N+1 \k) (= —EQ N \k)) close to k — it /a for several choices of c and c' in the 
case of TV = 4 (a) and N = 5 (b). Here, the hopping at the 1st leg and at the N-th leg ( see Fig. 1 ) 
are modified ascxt and c' x t, respectively. In each figure, the dotted line express the fitting by using 

E { N+1) (k)oc\k-7T/a\ N . 



dependence of the magnetic moment at the edge mix °c U l /^ N ^ should be observed even in the 
present situation if the anomalous power-law dependence of the magnetic moments and of the critical 
temperature discussed above are originated from the dispersion relation close to E/t = 0. In Fig. 7, we 
show the magnetic moments at the 1A site of the AF state at T = for the N = 4 and N = 5 system 
as a function of U/t for several choices of c and d . Here, the power-law dependence ttjia oc JJ l ^ N ^ 
is observed even if the hopping integrals at the edges are modified. The result is the strong evidence 
that the anomalous power-law dependence found in the present work is originated from the power-law 
dependence of the energy dispersion close to E/t = 0. We note that the magnetic moments at the 
edges does not satisfy the simple relation my = —itinb and are obtained as my > {tumbI m tne 
asymmetric case with c < d though the total magnetic moment in the unit cell vanishes. Even in such 
a case, as long as U/t <S 1, Im^l is also proportional to JJ 1 ^ N ~ l "> as well as my. 

4. Summary 

In the present work, we applied the mean-field approximation to the Hubbard model on the zigzag 
GNR and studied properties of the magnetic ordered states. The spontaneous magnetic moments at the 
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Fig. 7. Spontaneous magnetic moment at the 1 A site (upper zigzag edge in Fig. 1) m\A of the AF state with 
N — 4 and AT = 5atT = 0asa function of U/t for several choices of c and c'. The solid lines express 
the fitting by m 1A oc C/ 1 /^" 1 ). 

zigzag edges and the transition temperature of the AF state were investigated in detail as a function of 
on-site repulsion U for U/t <C 1. 

We can obtain the two kinds of the ordered states; one is the AF state satisfying = 
— m(jv-j+i)B and the other is the F state with rriiA = mnv-?+i)B- The AF state is more stable 
than the F state though the energy difference per one carbon atom between the two magnetic ordered 
states becomes smaller with increasing the width N. Therefore, the AF state was investigated in detail. 
Due to existence of the states localized around the zigzag edges close to Fermi energy, the magnetic 
moments at the zigzag edges are far bigger than the others for U/t <C 1, and show characteristic 
U/t dependence as mw oc (U ' /t) l ^ N ^ l \ Also, the transition temperature T Ci af shows the power- 
law dependence T c af oc {U /t) N /^ N ~ l \ which is analytically demonstrated from the divergence of 
the corresponding susceptibility. Discrepancy in the power of vti\a and T c af can be well understood 
by considering that the AF transition is dominated by the edge states and by assuming the simple 
dimension analysis, T Cj af oc U x txi\a- Therefore, we can conclude that the both anomalous U/t de- 
pendences are originated from the power-law divergence of the DOS close to Fermi energy. Actually, 
the power-law dependence of the magnetic moment at the edges are observed if the hopping integrals 
at the zigzag edges are modified, where the power of the DOS is unchanged. 
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